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ABSTRACT 



The amplitudes for decays of the type B^^s Dg ^Ds^-, have no factorizable contributions, 
O . while B^ s Dg^dD*^, and B^^g D*^^D*^^ have relatively small factorizable contributions 
through the annihilation mechanism. The dominant contributions to the decay amplitudes 
arise from chiral loop contributions and tree level amplitudes which can be obtained in terms 
p i' of soft gluon emissions forming a gluon condensate. We predict that the branching ratios for 
the processes 5^ ^ D+D;, 5^ ^ Df*D_- and B^ DfD;* are all of order (2 - 3) x 10-^ 
O : while 5° ^ DjD^, 5° ^ D+*Dj and S° ^ DjD^* are of order (4 - 7) x lO'^. We obtain 
branching ratios for two D*'s in the final state of order two times bigger. 



1 Introduction 



Recently, many theoretical and experimental studies in 5-meson physics have been done. The 
most extensive investigations have been done for cases of S-decay modes into the light mesons 
(e.g. B TTTT, B —>■ Kn) due to the their importance in the determination of Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. On the other hand 5-decays into charmed 
mesons also present an important issue for experimental and theoretical analysis. As it was 
pointed out by the authors of Pll2j, final state interactions (FSI) including charmed interme- 
diate states can give significant contributions to the decay amplitudes, especially in the case of 
the B — > Ktt decay modes In addition, branchig ratios for B^ decay into D^D^, D^*Df, 
£)-£)*+ and D*~D*~^ [2] have been measured. A previous theoretical study of the decay mode 
Bd — »• D~Df |lj showed that the non-factorizable contributions coming from chiral loops give 
an increase of order 10% in the rate compared to the factorized limit. The rate found in ^ 
is in good agreement with the experimental result. On the other hand, the B^ DfDj and 
B^ D^D^ decay modes have no factorizable amplitudes and they are realized only through 
non-factorizable contributions as it was shown in At the quark level, these decays a priori 
proceed through the annihilation mechanism bs — cc and bd cc, respectively. Within the 
factorized limit this mechanism will give a zero amplitude due to current conservation, as in 
the case of K^W [6:. Therefore 5,,, ^ D.^^D^^^ , 5^,, ^ D.^^Dl^ , and 5^,, ^ Dl^Dl^ 

present a fertile ground for the investigation of non-factorizable 1/Nc suppressed contributions. 

In this paper we want to extend our previous study of B^^g —>■ Ds ,iDs^d decays to cases 
with one or two D*'s in the final state. Namely, at 5-factories the decays to one pseudoscalar 
and one vector Z)-meson are easier accessable due to better statistical accuracy (they can 
be reconstructed more inclusively than decays to two pseudoscalars) jTj. In this case non- 
factorizable contributions arise due to the annihilation mechanism. However, its amplitude is 
suppressed by a numerically less favorable combination of Wilson coefficients and is expected to 
be of the same order of magnitude as non-factorizable contributions. Since the energy release 
for i?-decays to two charm mesons is relatively small (of order IGeV), the QCD improved 
approaches used to describe 5-decays to light mesons are not expected to hold. Therefore, 
we develop a different approach. 

As in [5^ our framework will be threefold: We use the standard effective Lagrangian approach 
for the quark process bq — > cc (where q = d,s) found at a scale below rric [9^ . Second, we use 
heavy-light chiral perturbation theory for interactions between heavy mesons and light pseudo- 
scalar mesons ^U] to calculate non-factorizable contributions in terms of chiral loops. Third, 
to estimate the contributions from 1/Nc suppressed terms at tree level |TT1 ^] within heavy- 
light chiral perturbation theory, we use a recently developed Heavy Light Chiral Quark Model 
(HLxQM) [in] based on the Heavy Quark Effective Field Theory (HQEFT) [Hj. 

Within our approach, the chiral symmetry breaking scale A^ ~ 1 GeV is the matching scale 
for pertubative QCD, chiral perturbation theory, and HLxQM. The latter is also the bridge 
between them. In the next section 2, we describe our framework. In section 3 we present the 
factorizable (annihilation), gluon condensate and chiral loop contributions to the amplitudes. 
In section 4 we give the numerical results and conclusions. 
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2 Framework 



2.1 Effective Lagrangian at quark level 

Based on the electroweak and quantum chromodynamical interactions, one constructs an effec- 
tive Lagrangian at the quark level in a standard well known way: 

i 

where Cj(/i) = ~^VcbV^ai{fi), q = d, s and aj(yu) are dimensionless Wilson coefficients that 
carry all information of the short distance (SD) physics above the renormalization scale /x. The 
matrix elements of the operators on the other hand, take care of all non-perturbative, 

long distance (LD) physics below /i. The relevant operators in our case are 

Ql = ^{IlI'^l) {cLlaCL) , Q2 = ^{cLl'^hi) {iLlaCL) , (2) 

where L denotes a left-handed particle. Contributions from other (say, penguin) operators are 
neglected due to smallness of their Wilson coefficients. 

In order to obtain all matrix elements of the Lagrangian we need the Fierz transformed 
version of the operators in ((21). To find these, we use the relation: 

^ij^ln = Jj-SinSlj + 2 tfj , (3) 

where i, j, k and n are color indices running from 1 to 3 and a is a color octet index. One 
obtains 

Ql = -^Q2 + 2Q2 , = ^Qi + 2Qi , (4) 

where the superscript F stands for "Fierzed", and 

Ql = 4(gL7"^"feL) {cLlaf^CL) , Q2 = 4 {cLl''t%L) (qL^afcL) , (5) 

where t"" denotes the color matrices. 

In order to calculate the matrix elements of the operators Qi and Q2 we will use a version 
of the Heavy Light Chiral Quark Model (HL^QM) developed in 13J. It belongs to a class 
of models extensively studied in the literature ^Hj - [2^ and is appropriate for describing 
interactions in which the transferred enegry is of the order 1 GeV. This sets the scale in to 
yLi ~ ~ 1 GeV, which is by construction the matching scale within our approach. At this 
scale one finds Oi ~ —0.35 — 0.07i and 02 — 1.29 -|- O.OSi IHIH]. Note that the a^'s are complex 
below the charm scale. 

2.2 Heavy light chiral perturbation theory 

The construction of heavy light chiral perturbation theory is based on the heavy quark effective 
theory (HQEFT) which is a systematic l/mq expansion in the heavy quark mass mg. The 
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Lagrangian is obtained by replacing the heavy quark Dirac field Q{x) = b{x),c{x) or c with 
a "reduced" field Q^^\x) for a heavy quark, and QIj^\x) for a heavy anti-quark. These are 
related to the full field Q{x) in the following way: 



(6) 



where P± = (1 ± 7 ■ f )/2 are projecting operators and v is the velocity of the heavy quark. The 
Lagrangian for heavy quarks then reads: 



(7) 



where is the covariant derivative containing the gluon field and 0{m'^) stands for the l/mQ 
corrections |tl2J, which will not be considered in this paper. 

After integrating out the heavy and light quarks, the effective Lagrangian for heavy and 
light mesons up to 0{mQ ) can be written as [T3t I22j: 



+ .... 



where -ff^^^ is the heavy meson field containing a spin zero and a spin one boson: 



(9) 



The field P^^\pj^^) (M = jj, for a vector and M = 5 for a pseudo-scalar) annihilates (creates) 
a heavy meson containing a heavy quark (anti-quark) with velocity v. Furthermore, = 
5baD^ — V^, where a, h are fiavor indices. The vector and axial vector fields and are 
defined as: 



(10) 



where / is the bare pion coupling, and n is a 3 by 3 matrix which contains the Goldstone 
bosons TT, 7] in the standard way. The ellipses in (jS)) denote terms of higher order in the 
chiral expansion. 

Based on the symmetry of HQEFT, we can obtain the bosonized currents. For a decay of 
the hq system we have fSl 122] : 



fill 



where an = fny/^H- In the limit itiq —>■ 00, = cxd = 01 h, but there are l/mq and 
perturbative QCD corrections to this limit [13 E]. Here Q^'^ is the heavy 5-quark field, is 

its velocity, and H^^^ is the corresponding heavy meson field. 
For the M/^-boson materializing to a Z)-meson we obtain: 



(-) 



an 



Tr 



(12) 
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Figure 1: Factorized contribution for 
not exist for decay mode B^ DfD~. 



D^Dg through the spectator mechanism, which does 



where ■* and v (= fg) is the heavy quark field and the velocity of the c quark, respectively, 
and ifg ^ is the corresponding field for the D-meson. Similar, for the b ^ c transition, the 
bosonized current is: 



-C(cu)Tr 



(13) 



where Q^^-* and Vc is the heavy quark field and the velocity of the c-quark, respectively. Fur- 
thermore, C(<^) is the Isgur-Wise function for the B ^ D transition, H^,'^'^ the heavy D-meson 
field(s), and u = ■ Vc = vi, ■ v = 1/(2k), where k = M^/Mb. 

For the weak current for DD production (corresponding to the factorizable annihilation 
mechanism) we obtain: 



(14) 



where X = v ■ Vc = (1/(2k^) — 1), and C(~'^) is a (complex) function less studied and not so 
well known as C(^) iii f)13|) . 

The factorized contributions for the spectator and annihilation diagrams are shown in Figs. 
Hand 121 The first diagram does not give any (direct) contributions to the class of processes we 
consider, but is still important because it is the basis for our chiral loops. 



The non-factorizable chiral loop contributions to the amplitudes can be visualized in Fig. 
where the heavy pseudo-scalar mesons are represented with single lines, the heavy vector 
mesons with double lines, the light kaons with dashed lines, and the weak vertex with two 
circles. As seen from Fig. El the chiral loop contribution can be divided into two topologies 
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Figure 2: Factorized contribution for — > DfD^ through the annihilation mechanism, which give 
zero contributions if both Df and D~ are pseudoscalars. 



that we will denote as topology I and topology II, respectively In the topology I (diagrams 
^01) ^11 ! ^21, ^22, ^31(^41)5 and ^32(^42) in Fig. ^ the _B-meson radiates K'^ becoming a bs 
state that decays into D~{D~*) and D'^{Df*). The meson is then reabsorbed by D~{D~*) 
giving the D'^D-*) state. In the topology II (diagrams A02, ^12, ^13, ^23, ^24, ^33(^43) , 
^34(^144), and ^35(^45) in Fig. El), the 5-meson decays into D-{D-*) and D+{D+*) which 
then re-scatter by means of a kaon into Dj{Dj*) and D'^{Df*). In both cases, the amplitudes 
for intermediate 5° D^D~ or Bg D^D^ processes (and similar for one or two vectors 
in a final state) can, within the factorized limit, be written as the product of currents ()12j) 
and (fT^ multiplied by the numerically favorable Wilson coefficient combination (C2 + Ci/Nc). 
Vertices describing absorption or radiation of a kaon are given by (jH)). There are of course 
also factorizable loop contributions, but these are included in the decay constants and the 
Isgur-Wise functions. 

The calculation of the chiral loop amplitudes includes the calculation of a divergent integral 
of a form: 

In the dimensional regularization method, the integral can be rewritten as: 



/(t^l,f2 



1-x' 



where g'^ is a D dimensional matrix tensor and 



Ii = 



d^k 



167r2 



A - In 



K 



+ 1 



(16) 



(17) 



{27t)d F - ml 

with A = 2 / {A — D) — 'Je + ln47r. Here r = r{x) with x = f 1 ■ f 2 is a function defined by: 

r(x) = ^ ^ ^ In (^\x\ + Vx'^ — 1^ , (18) 



for a; > 0, and 



for a; < 0. 



X 



==(ln + x'^ — 1 



(19) 
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Figure 3: Feynman diagrams for the chiral loop contributions. For S° Df*Dg * , the dia- 
grams A^i stand for contributions proportional to the Levi-Civita term. 
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2.3 The heavy hght chiral quark model (HLxQM) 

The HLxQM Lagrangian can be written as 

-^HLxQM = ^HQEFT + ^xQJ^I + ^Int ■ (20) 

The first term describes the interaction of heavy quarks in ((7j). The second term describes the 
interactions of hght quarks with hght (Goldstone) mesons in terms of the chiral quark model 

(xQM) [13 [131211231211: 

i^xQM = x[l^{^D^ + V^ + l5A^)-m]x. (21) 

Here m = 0.23 ± 0.02 GeV is the SU{3) invariant constituent light quark mass, and x is the 
flavor rotated quark fields given by xl = ^^Q'l and xr = ^Qr^ where = {u,d,s) are the 
light quark fields. The left- and right-handed projections ql and are transforming after 
SU{3)l and SU{3)r respectively. In we have discarded terms involving the light current 
quark mass which are irrelevant in the present paper (but become important when calculating 
counterterms) . The covariant derivative -D^ in ()21|) contains the soft gluon field forming the 
gluon condensates. The effects based on can be calculated by Feynman diagram techniques 
as in [HI 1121 Uni 1211 12n] or by the means of heat kernel techniques as in [THl 1211 121] • 

The interaction between heavy meson fields and quarks is described by [IHl dZl UHl 1201 HI] : 



^Int — —Gh 



(22) 



with the coupling constant Gh = y/'^rnp/ f , where p is a hadronic parameter (of the order one) 
depending on m [T^ . 

Within the model, one finds the following expression for the Isgur-Wise function jl3j : 

C(uj) = -^{l-p)+pr{uj), (23) 
1 + ui 

where r{u) is given by (fTH|) . In the simple expression (j^^ chiral loop and perturbative QCD 
corrections down to the scale p = A^, have not been taken into account. Nevertheless, it gives 
a good description of the Isgur-Wise function ^3]. For negative values of the argument the 
Isgur-Wise function C(~'^) is, within the HLxQM ^S], still given by (j23|) with uj —A. 

Performing the bosonization of the HLxQM, one encounters divergent loop integrals which 
will in general be quadraticly, linearly and logarithmicly divergent piHJ. The quadraticly and 
logarithmicly divergent integrals are related to the quark condensate and the bare pion decay 
constant / [2S] , respectively. The linearly divergent integral (which is finite within dimensional 
regularization) is related to the axial coupling (7^ in (|H]). 

The gluon condensate contribution to the amplitudes can be written, within the framework 
presented in the previous section, in a quasi-factorized way as a product of matrix elements of 
colored currents: 



{DtD;\Cwm%F = 8C2{DtD;\cL7^t''cL\G){G\dLYt''bL\B'^\ 



') , (24) 
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Figure 4: Non-factorizable contribution for — >■ D^D~ through the annihilation mechanism 
with additional soft gluon emision. The wavy lines represent soft gluons ending in vacuum to 
make gluon condensates. 



where G in the bra-kets symbolizes emission of a gluon as visualized in Fig. EJ The left part in 
Fig. 0] gives us the bosonized colored current: 



^6 Jig 



GuuTr 
647r 



ia^'' - F {a^\ ^ ■ v,}) , (25) 



where G^^, is the octet gluon tensor, and F = 2T^f'^/{m? is a dimensionless quantity of 



the order 1/3. The symbol { , } denotes the anti-commutator. 

For the creation of a DD pair in the right part of Fig. 0] (the analogue of (j25|) ) one gets: 



IG 



G"^ Tr 

1287rm(A - 1) 



X {Xa^"^ + {a'^^7■ At;} 



where 



X^-(A-l)r(-A) 



TT 



(26) 



(27) 



and Av = Vc — v . Multiplying the currents in Eqs. and ^I^i. and using the replacement: 



,a 



9sG°^uGa(3 ^ 47r ( G ) — igi,a9ul3 - 9tif39iya) , 



(2J 



we obtain a bosonized effective Lagrangian term which is l/N^ suppressed compared to the 
factorized contributions. This effective Lagrangian term corresponds to a certain linear com- 
bination of a priori possible l/N^ suppressed terms at the tree level in the chiral perturbation 
theory sense. 



3 Calculating the decay amplitudes 

First, for comparison, we give the factorized amplitude for the process 5° D^D's 



A{B<^ ^ D+D;)f = {G2 + -f) C(^)/dMb VMsMz) {\ + uj 



N, 



(29) 



The MM, M = D'^* decay amplitudes can be written in the following form: 

A{B DD) = F^^, A{B D*D) = iF^^ sd* ■ v , (30) 

AiB DD*) = iF^^ Ed, ■ , (31) 

A{B ^ D*D*) = Fo^^ Ed, ■ e^. + F^"" ed- ■ v eq, ■ + F^ 2e„^^,<t;^£^.£^-, , (32) 

with the reduced amplitudes Ff'^ , containing the contributions coming from the tree level 
factorizable contributions -Fjj^, gluon condensates Fj^^^ and chiral loops Ff^J^: 

where N = P,V mean pseudoscalar D and vector meson D*, respectively. 

The factorizable amplitude for B^ MM, M = _D+,_D+* comes from the annihilation 
diagram: 

{D;Di\Cw\W)p = 4(Ci + ^C,){D;Dt\cEi,CL\0){0\d2YbLm , (34) 

and is proportional to a numerically non-favorable combination of Wilson coefficients. Using 
Eqs. (fTT|) and (fT^ we obtain the following values for the reduced amplitudes: 



Fojc = -Flfc = F^Jc\l^^ = 2(Ci + ^)C(-A)a^«: ^MsMMr , (35) 

FqJc = FqJc = Fijc = • (36) 

Using (25) - (28), we obtain the gluon condensate contributions (the color suppressed l/Nc 
contributions at the tree level) to the reduced amplitudes illustrated in Fig. HI 



F,^l = 3S (x + ^iX-1)) , Fl^, = S {X[l-2F]+A[\ + 2F]) , (37) 



= ^ (X [3 + 2F] + 4 [(A - 2) - 2F]) , F^ = 2k'SX (A + 1) , (38) 

F^^, = -2^S (X - 4) , F^l = 2k' S {X - 4)(1 + 2F) , (39) 
where F is defined below (jH)), X is defined in (E7j) . and 



C2 {GhVME 



3 



3-29m(A-l) n ' ^ ' 

In the evaluation of the integrals (fTH|) and p7|l appearing in the reduced chiral loops ampli- 
tudes Fflf , we use the MS scheme and we take yU ^ ~ 1 GeV. The integral (fTTj) contains 
a logarithmicly divergent term and a constant term. However, additional contributions to the 
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constant term might come from counterterms. The finite part of these terms is unknown and 
therefore, in our numerical computation, we take into account logarithmic terms only, which are 
independent of the counterterms contributions, and we consider constant term contributions as 
theoretical uncertainities of our approach. 

However, the products of two Levi-Civita terms enter in our computation of the amplitudes 
coming from diagrams , , ^32 ? ^35 and ^445. Whitin dimensional regularization, these 
products are not uniquelly defined. This is related to the problem known in the literature as 75 
scheme dependance |j27j . However, this scheme dependence appears in the constant terms only. 
To estimate its infiuence on numerical results, we use two different appraches to the products 
of two Levi-Civita symbols. (One is the dimensoinal reduction |2Zj and the second one is a 
variation of the dimensional reduction in which the products of metric tensors in D dimension 
and the metric tensor in 4-dimensions are fixed on the 4-dimensional space.) 

The logarithmic contributions to the reduced amplitudes are: 



where 



Kci = ^K^^^[-2{l + uJ){r{-uJ)+r{-\)) + 2{uJ + \)], (41) 
Fl!, = -rK[~2K{r{-uj) + r{-X)) + 2{K + l)], (42) 



F^X = -^K[2K{r{-uJ)+r{~X)) + 2{K + l)], (43) 
FlI = ^K^J^{2{u; + l)), (44) 



FU = zK^^{{K + l)2[r{-uj)+r{-X)]~2{u; + X){H + G)~2K'), (45) 
^ -^/f^., (46) 



K ^ \C, a^)an (^-f\ ' II ^M^Mn* Mj, , h^^\n^, (47) 



(_ri-u)+ul (l + urj-u)) (1 + Ar(-A)) 

G = K K , h = — . (48) 

1 — uj'' 1 — 1 — A"^ 

Note that is a priori proportional to the Wilson coefficient combination (C2 + Ci/Nc) , as the 
amplitude in (j^ . But because the factor 1//^ is already of order l/NcWe replace {C2+Ci/Nc) 
by just C2 in (gTj). 



4 Results and discussion 

In our calculation we used the following input parameters: as = old = 0.33 GeV^'^^^, p = 1.05, 
Gh = 7.5 GeV-1/2 and {^G^} = [(0.315 ± 0.020) GeV]^ 9a = 0.6 [23, U = 0.093 

GeV and K = 0.37. 
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i 


X 10^ 


Fir X 10^ 


FlI X 10 ^ 


Fir X 10^ 


Frr X 10^ 


Fl/ X 10^ 


fc 





0.12i 


-0.12i 








-0.12i 


gc 


-0.17-0.71i 


1.11 + 0.14i 


-0.51 + 0.87i 


0.13-0.22i 


-0.13+0.06i 


0.22 -O.OQi 


d 


0.91-1.20i 


-0.23 + 0.19i 


-0.08 - 0.21i 


0.56+0.03i 


-0.24-0.45i 


-0.09 



Table 1: Reduced amplitudes for the S° MM, M — Df, D+* decay modes. The results are 
given in GeV. 



i 


Fir X 10-^ 


Fir X 10-^ 


Fir X 10-^ 


Fir X 10^ 


X 10^ 


Fir X 10^ 


fc 





0.511 


-0.511 








-0.511 


gc 


-0.74-3.16i 


4.95 +0.631 


-2.27 +3.691 


0.60-1.001 


-0.59+0.251 


0.99 -0.421 


cl 


3.78-4.961 


-0.96 +0.781 


-0.34 - 0.851 


2.13+0.101 


-0.91-1.721 


-0.34 -0.021 



Table 2: Reduced amplitudes for the —>■ NN, N = D~^,D~^* decay modes. The results are 
given In GeV. 



The reduced amphtudes for S° MM, M — Df, D^* are presented In Table 1. We find 
the following branching ratios: 

D+D-) 





- Dtor 


= 2.5 X 10-^ , 


5r(5° 


BriB"" - 




= 3.3 X 10-^ 


Bt{BI 


BriB"" - 




= 2.0 X 10-^ , 


Br{Bl 


Br{B° 




= 5.4 X 10~\ 


Br(B', - 



4.5 X 10-^ , 


(49) 


6.8 X 10-^ , 


(50) 


4.3 X 10-^ 


(51) 


9.1 X 10~^ 


(52) 



D*+D~*) 

The contribution of the constant term and the corresponding countcrterm can change the 
branchig ratio for i?-meson decaying into two pseudoscalars by about 10%, while in the case 
of decay into one pseudoscalar and one vector D-meson, this contribution is in the range of 
20 — 40%. In the case of S-meson decaying into two vector mesons, the constant term is 
estimated to be 2-8 times larger than the logarltmic contribution, depending on the choice of 
the scheme in which the products of two Levi-Civita terms are considered. The uncertainty in 
input parameters can result in an additional error for the branching ratios. We estimate that 
it can be of the order of 20%. Within our approach the l/mg corrections, with Q — c,b have 
been neglected. At least the l/rric corrections might be important. 

The study of dominant contributions in the S° F)fD~, 5° F)~^*F)~ and F)fD~* 
decay amplitudes is very important for our undcrstaning of the color suppressed contributions 
to the 5-mesons decaying to two charm mesons. This is even more important knowing that 
the experimental rates for B^ decay into D~Dr D~*Dr D'D*'^ and D*~D*'^ have very 
small color suppressed contributions and therefore the decay amplitudes we consider open a 
window for studies of color suppressed effects in 5-decays to two charm mesons. The chance 
for experimental measurements of these decay rates at i?-factories makes the study of the decay 
mechanisms in these decays even more important. 
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5 Appendix A: Chiral loops with "superpropagator" 

In the following, we will shortly present a "superpropagator method" which enables us to 
calculate all contributions of one kind of topology to all decays of a type MM, M = 

D*,Df* as one compact calculation. This means that instead of calculating all the diagrams 
in Fig. 01 we only need to calculate two contributions, each one coming from one topology. 

The superpropagator is a propagator of the H^'^^x) field and therefore includes propagators 
of both heavy pseudo-scalar and heavy vector meson. It is defined in a standard way as a 
contraction of H^^\x)ap and H'^^\y)i^\, where a, /3, k, A are Dirac spinor indices. In momentum 
space, we obtain the superpropagator 



c{±) 



{k) 



1 



2{±v -k + it) 



where k is the momentum and 



(53) 



(54) 



For the effective Lagrangian of topology I we have: 
Ci{B DD) = 



13a 



X 



while the effective Lagrangian of topology II can be written as: 



CniB ^ DD) 



X 



(+) 



(+) 



I3a 



af3;KX 



(+) 



Xk 



From and (j37|) one can derive the (reduced) amplitudes already given in pTjl - ljlUj) . 



(55) 



(56) 



(57) 
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